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ABSTRACT 

In the current era of high-precision CMB experiments, the imprint of gravitational lensing 
on the CMB temperature is exploited as a source of valuable information. Especially the 
reconstruction of the lensing potential power spectrum is of great interest. The reconstruction 
from the optimal quadratic estimator of the lensing potential, though, is biased. As long as 
the intrinsic CMB fluctuations are Gaussian this bias is well understood and controlled. In 
the presence of non-Gaussian primordial curvature perturbations, however, the CMB also 
acquires a non-Gaussian structure mimicking the lensing signal. Concentrating on primordial 
non-Gaussianity of local type, we address the resulting bias by extracting the lensing potential 
power spectrum from large samples of simulated lensed CMB temperature maps comprising 
different values of /nl- We find that the bias is considerably larger than previous analytical 
calculations suggested. For current values of /nl and a sensitivity like that of the Planck 
mission, however, the bias is completely negligible on all but the largest angular scales. 

Key words: cosmology: cosmic microwave background, gravitational lensing, methods: nu- 
merical 



1 INTRODUCTION 

The cosmic microwave background (CMB) carries an incompara- 
ble content of cosmological information. After removing all fore- 
grounds its tiny temperature fluctuations observed today mirrors the 
physical conditions in the Universe approximately 400,000 years 
after the big bang. The primary tool to access this information is the 
CMB temperature power spectrum. Additionally, CMB polariza- 
tion is an important diagnostic of gravitational waves. Gravitational 
lensing of the CMB by the intervening large scale structure, how- 
ever, alters the shape of the observed power spectra (see Lewis & 
Challinor 2006, for a review). In particular, it inevitably introduces 
B-modes confusing the imprint of primordial tensor perturbations. 
Therefore, in the context of parameter analysis and gravitational 
wave detection CMB lensing can be regarded as contaminant. 

However, the unprecedented sensitivity of ongoing and future 
CMB experiments allows for a new perspective: The lensing sig- 
nal in the CMB temperature itself contains valuable cosmological 
information. Precise knowledge of the lensed CMB power spec- 
trum can be used to improve on parameter estimations from the 
primary CMB breaking several degeneracies (Hu 2002; Stompor & 
Efstathiou 1999; Smith et al. 2006). Furthermore, lensing probes 
the (integrated) matter distribution from today's observer up to the 
last scattering surface. As the growth of large scale structure is sen- 
sitive to neutrino mass so is the lensing signal thereby tightening 
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up the mass limits (Lesgourgues et al. 2006). Besides these appli- 
cations in cosmological parameter studies, direct reconstruction of 
the lensing effect is also of major importance. In order to access the 
primordial B-mode signal the lensing effect needs to be removed 
from the observed CMB requiring the lensing potential (Knox & 
Song 2002; Seljak & Hirata 2004). 

There are several approaches to extract the lensing potential 
from the observed CMB (Zaldarriaga & Seljak 1999; Guzik et al. 
2000; Hu 2001b; Hirata & Seljak 2003; Okamoto & Hu 2003). 
Okamoto & Hu (2003) is very well suited for the instrumental 
noise level of current experiments regarding precision and numeri- 
cal complexity. Being quadratic in the lensed temperature field this 
estimator exploits the lensing induced non-diagonal elements of the 
CMB covariance matrix. In the next step the reconstructed lensing 
potential can be used to estimate its power spectrum. The resulting 
estimator, however, is biased. There are several contributions to this 
bias, which can be classified according to the order of the lensing 
potential involved. The zeroth order contribution is dominant and 
given by the Gaussian variance of the estimator of the potential. 
Smaller but substantial contributions are of first and second order 
in the lensing potential and arise from the non-trivial part of the 
(lensing induced) CMB trispectrum. Up to date all relevant contri- 
butions to this so-called reconstruction bias are identified and under 
control (Kesden et al. 2003; Hanson et al. 201 1). 

Since not only the information about the lensing potential 
power spectrum itself, but in part also the reconstruction bias, orig- 
inates from the connected part of the lensing induced CMB four- 
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point function, any kurtosis present in the observed CMB will af- 
fect the power spectrum extraction, too. A non-trivial four-point 
function may arise from primordial non-Gaussian curvature per- 
turbations, which are predicted by all inflationary models. Current 
observations suggest that the deviations from Gaussianity are quite 
small. Nevertheless, depending on the inflationary scenario, they 
can be rather sizable. 

In this paper we extend the analytical work of Lesgourgues 
et al. (2005) to the full sky and investigate the impact of an in- 
trinsic non-trivial CMB trispectrum on the lensing power spectrum 
reconstruction using numerical simulations. 

The structure of this paper is the following: in Section 2 
we briefly review the key aspects of CMB lensing power extrac- 
tion. Then, in Section 3, primordial non-Gaussian curvature per- 
turbations of local type and their connection to a non-trivial CMB 
trispectrum are explained in detail. Our method to determine the 
additional reconstruction bias in the presence of primordial non- 
Gaussianity is presented along with the results in Section 4. Finally, 
Section 5 is devoted to a short summary of the main findings of this 
work. 

As reference cosmology we choose the parameters compiled 
in the WMAP5+BAO+SN data set 1 . The matter content is de- 
scribed by Q. c h 2 = 0.1143 and Q b h 2 = 0.02256, while the cosmo- 
logical constant is given by Q A = 0.721. Primordial perturbations 
are characterized by the scalar amplitude A = 2.457 • 10~ 9 and spec- 
tral index n s = 0.96. The value of the Hubble constant measured 
today is H = 100/ikms" 1 Mpc~' with h = 0.701 and the optical 
depth is assumed to be r = 0.084. 



where 



E c ,..„ = ^a(a + 1) • ■ ■ n(n + 1), n fl ...„ = V(2«+ 1) - - - (2/z + 1). (5) 

The sixfold indexed quantities in braces denote Wigner 3 j symbols 
(see Abramowitz & Stegun 1972). 

Okamoto & Hu (2003) realized that it is possible to extract 
the lensing potential from equation (3) by taking the correlator of 
the lensed multipoles while keeping the realization of the lenses 
fixed. This procedure returns the multipoles of the lensing poten- 
tial weighted by a function just depending on the intrinsic CMB 
power spectrum. Since lensing correlates temperature modes across 
a whole multipole band (whose width is given by the corresponding 
band power in the deflection field) one can construct an estimator 
of the lensing potential as weighted sum over lensed multipoles: 
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Here we have slightly changed our notation. From now on mul- 
tipoles with a tilde denote the observed temperature fluctuations 
containing the lensing signal, as well as all kinds of experimental 
noise. The weighted sum needs an appropriate normalization given 
by 
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Under the assumption that both the intrinsic temperature fluctu- 
ations and the lensing potential are homogeneous and isotropic 
Gaussian random fields, i.e. 



2 CMB LENSING POWER SPECTRUM 
RECONSTRUCTION 

On their way from the last scattering surface to today's observer 
CMB photons are deflected by the intervening large scale struc- 
ture. This gravitational lensing effect on the CMB is well de- 
scribed by a simple remapping of the temperature contrast, 0(n) = 
Ar(n)/r CM B, 



0(n) = [n + V0(n)] 



(1) 



where a tilde denotes lensed quantities. The transition from the 
original direction on the sky to the lensed one is mediated by the an- 
gular gradient of the lensing potential <p(n), i.e. the projected New- 
tonian gravitational potential ^(n) 
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where optimal means that they minimize the estimator's Gaussian 
variance. The spectrum Cf®. denotes the observed one, i.e. it con- 
tains the lensing signal and all kinds of instrumental noise and fore- 
grounds. 

Given the estimated multipoles of the lensing potential field 
one can proceed to estimate the corresponding power spectrum via 
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This integral is formulated for a flat universe where the angular 
diameter distance coincides with the comoving distance^-. It ranges 
to the last scattering surface at comoving distance %* . 

In harmonic space and linear in the deflection field the effect 
of lensing can be expressed as (Hu 2000) 
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The parameter values can be obtained from 
http : //lambda . gsf c . nasa . gov/product/map/dr 3/parameters . cfm. 



Accordingly, the expectation value of the reconstructed lensing 
potential power spectrum involves the four-point function of the 
lensed temperature contrast 
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This estimator, however, is biased. While the connected part of the 
four-point function contains the information about the true lensing 
potential power spectrum, its Gaussian part gives the largest contri- 
bution to the bias which is identical with the normalization of the 



estimator, i.e. N' 
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A ; . In addition, there are two more terms that 
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nced to be considered. The first one is linear in the lensing poten- 
tial spectrum. It was first derived by Kesden et al. (2003). Since it 
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Figure 1. Biases involved in the reconstruction of the CMB lensing po- 
tential power spectrum ', starting with the optimal quadratic estimator 
of the lensing potential defined in equation (6). The terms N) ' and 
are analytical estimates of the additional bias arising from primordial non- 
Gaussian curvature perturbations of local type where /nl = 100 has been 
assumed. 



is only substantial on small scales a treatment in the limit of a flat 
sky is well applicable 
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with the additional constraint li + It = L = I3 + I4. The explicit 
expressions for the flat-sky analogues of the normalization and 
weights are given in Appendix A. 

Furthermore on large scales Hanson et al. (2011) identified a 
negative bias which is quadratic in the lensing potential spectrum. 
Its approximative formula 
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is in excellent agreement with results from numerical simulations. 
All three biases are shown in Figure 1 . 



3 LOCAL NON-GAUSSIANITY 

In the derivation of the optimal reconstruction weights and the three 
different biases of the estimated lensing potential power spectrum 
the intrinsic CMB fluctuations were always assumed to be Gaus- 
sian distributed. In this case, and under idealized experimental con- 
ditions, the non-Gaussian structure of the observed CMB solely 
arises from the lensing effect. Non-Gaussianity of the lensed CMB 
then manifests itself in the connected parts of the n-point correla- 
tion functions. Odd moments are not generated by lensing as long 
as both fields are uncorrelated. Correlations, however, do arise on 
large scales due to the integrated Sachs-Wolfe (ISW) effect: the cos- 
mic large-scale structure does not only deflect the CMB photons 
but its late-time evolution also gives rise to an additional (large- 
scale) temperature anisotropy. On small scales there are correla- 



tions because of the non-linear ISW (or Rees-Sciama) effect, as 
well as the Sunyaev-Zel'dovich effects. These correlations result 
in a non- vanishing bispectrum of the lensed CMB. Thus, there are 
non-Gaussian signatures present in the lensed, i.e. observed, CMB; 
even when starting from Gaussian distributed temperature fluctua- 
tions. 

Inflationary models, however, predict at least small devia- 
tions from pure Gaussian initial conditions. While simple infla- 
tionary scenarios generate a nearly Gaussian spectrum of the pri- 
mordial curvature perturbations (Guth 1981; Bardeen et al. 1983; 
Mukhanov et al. 1992), more involved models are expected to in- 
troduce a significant amount of non-Gaussianity. These models can 
be classified by the shape of their bispectrum: besides equilateral 
(Alishahiha et al. 2004; Arkani-Hamed et al. 2004; Silverstein & 
Tong 2004; Seery & Lidsey 2005; Chen et al. 2007; Cheung et al. 
2008; Li et al. 2008) and orthogonal type (Senatore et al. 2010), 
local non-Gaussianity is an important category (Lyth et al. 2003; 
Dvali et al. 2004; Barnaby & Cline 2008; Byrnes et al. 2008; Al- 
abidi et al. 2010; Elliston et al. 201 1). 

Primordial non-Gaussian curvature perturbations of local type 
can be described starting from a Gaussian fluctuation field Ol(x) 
and setting 

* nl(x) = *l(x) + /nl (f L (x) - («£(*))) ■ (14) 

The amplitude of the non-Gaussianity observed is then set by the 
dimensionless parameter / NL which in this case is assumed to be 
constant although models with scale dependent amplitude do exist 
(Byrnes et al. 2010). 

The ansatz of equation (14) can be extended by adding a fur- 
ther term of the form ^nl < I ) l (x). We ignore this term in our analysis 
for simplicity and because g NL is observationally almost completely 
unconstrained. However, in order to be competitive with the /nl 
contributions, its value must be of the order y#NL ~ /nl but is not 
constrained by theory, in fact it is possible to tune the inflationary 
model arbitrarily in these parameters (Sasaki et al. 2006). Likewise, 
terms proportional to tnl, which are similar in magnitude relative 
to /nl because of the Suyama-Yamaguchi relation tnl > (6/5 /nl) 2 
(see, e.g. Smith et al. 201 1 ; Suyama & Yamaguchi 2008; Sugiyama 
et al. 2011; Suyama et al. 2010) are neglected. Neglecting these 
contributions, however, should at most result in a slight underes- 
timation of the impact of local primordial non-Gaussianity on the 
lensing power reconstruction. 

By construction, the four-point function of the curvature per- 
turbations acquires a non-trivial part. This is also true for the intrin- 
sic CMB fluctuations since they are simply related to those in the 
curvature field via an integration over the corresponding radiation 
transfer function A t (k) 
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with N L(k) being the Fourier transform of equation (14). 

The non-trivial part of the angular trispectrum of the CMB 
temperature contrast is most conveniently expressed in its fully re- 
duced form (at leading order in /nl) 
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(see Okamoto & Hu 2002, for details on this formalism). In the 
expressions above, /\m>(/c) is the power spectrum of the curvature 
perturbations and j ( (r) denotes the f-th spherical Bessel function 
of the first kind (see Abramowitz & Stegun 1972). 

The computation of the trispectrum given in equation (16) is 
numerically quite cumbersome. Okamoto & Hu (2002), however, 
found a very useful approximation formula by extending the Sachs- 
Wolfe approximation, i.e. 
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valid at small I into the acoustic regime 

' f 3 f 4 ( L ) ~ ibh hLh h huJm^L c e 2 c f 4 
with 

^ = J; f ^ P ™( k )jl [k(l0 - 7rec)] 

where rjo - »7rec is the conformal time spread between today and 
recombination. This approximation holds as long as the last scat- 
tering surface can be considered as thin and its temperature fluctu- 
ations only vary slowly. 



4 POWER RECONSTRUCTION BIAS 
4.1 Zeroth order calculation 

Since the first and second order bias of the reconstructed lensing 
potential power spectrum arise from the connected part of the (lens- 
ing induced) CMB trispectrum, one expects that also an intrinsic 
(non-trivial) trispectrum of the unlensed CMB contributes an ad- 
ditional bias. This bias was first addressed by Lesgourgues et al. 
(2005). In their derivation they assumed that both of the effects 
leading to a connected part of the observed CMB four-point func- 
tion are completely separable. Thus, they did not consider a so-to- 
speak lensed trispectrum but just added the connected part of the 
unlensed CMB trispectrum. In this approximation, and in the limit 
of a flat sky, the additional bias is given by 
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with the additional constraints li +I2 = L and I3 + I4 = L. The func- 
tional form of the flat-sky expressions are relegated to Appendix A. 
Lesgourgues et al. (2005) also computed the bias contribution from 
the g NL term which is not considered in this work. 

Starting from equation (11), the all-sky generalization of equa- 
tion (24) is 
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The quantity in curly brackets denotes a Wigner 6 j symbol (see 
Abramowitz & Stegun 1972). Using the approximative formula for 
the trispectrum of unlensed temperature fluctuations given in equa- 
tion (22) the similarity with the lensing induced CMB trispectrum 
(Hu 2001a) 
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is striking. There are two important consequences. First, the term 
involving the Wigner-6 j symbol, hereafter denoted by Nj^, should 
only contribute on small scales, so it can safely be treated in the flat 
sky-limit, i.e. by using the second part of equation (24) . Second, 
the other term, afterwards denoted by N { *\ should be proportional 
to C^ w implying a scaling of the form 
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for a scale-invariant power spectrum (cf. Liddle & Lyth 2000). 

In Figure 1 the two additional bias contributions A^ A) and 
N ( P are shown. The agreement with the results of Lesgourgues 
et al. (2005) is very good although they used the exact temper- 
ature transfer functions instead of the extended Sachs-Wolfe ap- 
proximation applied in this work. Only in the acoustic regime the 
N*~ L ' term shows substantial differences. There, however, the bias 
is completely dominated by the ' term, which is accurately de- 
scribed by the extended Sachs-Wolfe approximation even on small- 
est scales as the comparison with Figure 1 of Lesgourgues et al. 
(2005) shows. On large scales, the full-sky calculation does not 
greatly improve the results of Lesgourgues et al. (2005) obtained 
neglecting the curvature of the sky. For both methods one recov- 
ers the expected scaling relation with the multipole order (cf. again 
Figure 1 of Lesgourgues et al. 2005). 

As already pointed out by Lesgourgues et al. (2005) for cur- 
rent values of /nl, -10 < /nl < 74 (Komatsu et al. 201 1), the addi- 
tional bias is not a severe contaminant of the reconstructed power 
spectrum. On almost all scales the bias is more than four orders of 
magnitude smaller than the biases arising from the lensing induced 
trispectrum. Only on largest scales the contributions from the 
term can be substantial. Consequently, as long as a fully separate 
perturbative treatment of both, gravitational lensing, as well as pri- 
mordial non-Gaussianity, is applicable, the additional bias is well 
controlled and for current / NL parameters negligible. 

However, the validity of this separation ansatz is not obvious. 
In the presence of primordial non-Gaussianity the mode coupling 
due to lensing is more involved since here the perturbative expan- 
sions of the lensed temperature field in powers of the lensing poten- 
tial on the one hand and the unlensed temperature field in powers 
of the non-Gaussian curvature perturbations on the other hand com- 
pete. In order to circumvent the assumptions and limitations which 
are necessary for an analytical study of this problem we will resort 
to numerical simulations. Here all the different effects involved, pri- 
mordial non-Gaussianity, lensing, instrumental noise, can be incor- 
porated in a consistent manner. 
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4.2 Numerical simulations 

The starting point of the simulations is the large sample of non- 
Gaussian CMB realizations provided by Eisner & Wandelt (2009). 
The non-Gaussian structure arises from primordial curvature per- 
turbations of local type detailed in Section 3. The maps of Eisner & 
Wandelt (2009) are well suited for our analysis because they did not 
focus on the intrinsic CMB bispectrum but produced simulations 
with accurate non-Gaussian statistics to all correlation orders. Ev- 
ery realization consists of a set of Gaussian (©]?,„) and non-Gaussian 
(0^) multipoles, so that a CMB map containing the desired level 
of non-Gaussianity characterized by /nl can be built via 



®lm - ©L + /nl©^- 



(29) 



The cosmological background of these simulations is adapted from 
the WMAP5+BAO+SN data as summarized at the end of Section 1 . 

Having a non-Gaussian CMB map at hand, this can be lensed 
in the framework of LensPix 2 (Lewis 2005). The remapping tech- 
nique implemented in LensPix assumes the lensing potential being 
fully described by its power spectrum. Obviously, this is no longer 
true in the presence of primordial non-Gaussianity. However, for 
the reconstruction scheme detailed in Section 2 the CMB and the 
lensing potential are treated as statistically independent fields and 
higher order statistics of the lensing potential are not involved. 

Since the map is already given in multipole space we use a 
slightly modified version of LensPix which skips the simulation of 
a (Gaussian) CMB realization and starts directly from the temper- 
ature multipoles constructed according to equation (29). In order 
to get a lensed CMB map at HEALPix (Gorski et al. 2005) resolu- 
tion N Mc = 2048 which is cosmically accurate up to £ max = 2750 



one needs input multipoles up to fj, 



nput 



3000. The non-Gaussian 



maps of Eisner & Wandelt (2009), however, cover only a multi- 
pole range up to £f NL = 1024. We therefore extend the maps in 
the range ff NL < ( < 3000 with pure Gaussian multipoles. Thus, 
equation (29) needs to be replaced by 



0L+/nl©^, if £< 1024 
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The power spectrum needed for the construction of Gaussian CMB 
realizations has been computed using CAMB 3 (Lewis et al. 2000). 

This method is not expected to crucially influence the main 
conclusion of this work. The non-Gaussianity of the maps con- 
structed according to this prescription will be even milder due to the 
fact that scales smaller than about ten arcminutes are not affected 
by the primordial non-Gaussian curvature perturbations. Likewise, 
if there are serious concerns about this method, one can just think 
of the resulting maps as not being subjected to primordial non- 
Gaussianity of the local type but subjected to some intermediate 
type which is scale dependent. In any case, the influence of the re- 
sulting non-trivial intrinsic CMB trispectrum on the reconstructed 
lensing power needs to be studied. 

Finally, the lensed maps are exposed to instrumental noise 
only. But neither any further secondary anisotropies nor fore- 
grounds nor experimental complications like incomplete sky cover- 
age are included. Adapting the specification of the Planck satellite 
mission, the maps are first convolved with a Gaussian beam with 
°"fwhm = 1 arcmin. Before deconvolution, white Gaussian pixel 
noise with standard deviation cr N = 27 /jK arcmin is added. The 
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Figure 2. Illustration of the performance of the CMB lensing power re- 
construction discussed in the text and condensed in equation (35). The red 
boxes indicate the width of the multipole band over which the extracted 
power spectrum has been binned, as well as the ltr error for one realization. 
In addition, the reconstructions in presence of primordial non-Gaussianity 
for two different values of /nl are shown. Note that on large angular scales 
the bias for /nl = 1000 exceeds the range of the plot by far (cf. Figure 4) 
and is not shown for clarity. 



power spectra of the final maps are well described by the analytic 
expression derived by Knox (1995) 
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4.3 Efficient estimator 

For practical purposes, the estimator in the form of equation (6) 
is not suitable. One should use its angular space representation in- 
stead 
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Here the spin-gradient derivative V just acts on the spin part of the 
corresponding fields, the inverse weighted temperature map 



Cm ^f.obs 



(33) 



and the Wiener reconstruction of the unlensed temperature field 

(34) 
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[in "-f.obs 

(Okamoto & Hu 2003). This representation allows for the use of 
fast harmonic transform algorithms on the sphere. 

4.4 Results 

In order to determine the bias associated with the connected part 
of the trispectrum of the unlensed CMB which is generated by pri- 
mordial non-Gaussian curvature perturbations of local type we pro- 
ceeded as follows 

(1) construct N CMB maps following equation (30) setting 
/nl = 

(2) use these maps to build a sample of lensed maps which con- 
tain the desired instrumental noise properties as detailed in Sec- 
tion 4.2 
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(3) apply the efficient estimator of Section 4.3 to the sample and 
reconstructing the lensing potential power spectrum for each map 
according to 
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where the upper index runs from 1 to N. The three bias contribu- 
tions in the equation above have not been computed separately for 
each individual realization. Instead, the theoretical power spectrum 
of the lensed and unlensed CMB, as well as of the lensing poten- 
tial, again computed with CAMB, were used. To illustrate that in 
this way the (purely lensing induced) biases are very well captured, 
we plot in Figure 2 the lensing potential power reconstruction for 
/nl = obtained from hundred CMB realizations. The boxes in- 
dicate the width of the multipole band and the one sigma error for 
one realization. 

(4) compute the sample mean 



j_ V r <*<4,(0 



(36) 



(5) repeat steps (1) to (4) but now setting /nl = 100 in the first 
step 

(6) repeat step (5) but now starting with /nl = 1000 

(7) compute the bias for both values of /nl according to 
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It needs to be stressed that in steps (1) to (3) for all three values of 
/nl the simulation conditions have been the same. In other words, 
we kept track of the different random seeds involved in the map 
making procedure. Thus, the j-th map is always lensed by the same 
realization of the lensing potential and always contains the same 
realization of instrumental noise. This ensures that the differences 
in the reconstructed power spectra, i.e. the bias, are solely caused 
by the variations in /nl- The numerical analysis of this work has 
been earned out setting N = 100. 

In Figure 2 the reconstructed power spectra for both values of 
/nl are plotted, while Figure 3 and Figure 4 shows the correspond- 
ing biases. For comparison, also the bias derived by Lesgourgues 
et al. (2005) (cf. equation 24 and 25) is plotted. While the functional 
form of the bias found in the simulations resembles that of the ana- 
lytical zero-th order prediction, its amplitude is much larger. Inter- 
estingly, the disagreement in amplitude is different for the two parts 
of the analytic bias formula. For the smaller value of /nl the N 1 ^ 
part which gives the relevant contribution on large scales, needs to 
be scaled by a factor of roughly 100 to reach the same level as the 
bias found in the simulations. Its small scale counterpart, the A^ S) 
term differs approximately by a factor of 125. 

Looking at the bias obtained for /nl = 1000 shown in Fig- 
ure 4, one immediately realizes that in contrast to the analytical 
formula of Lesgourgues et al. (2005) the bias extracted from nu- 
merical simulations does not scale quadratically with /nl- This is 
particularly true for large angular scales. Here the agreement with 
the analytic expression, now phenomenologically amplified by a 
factor of 50, is much better than for /nl = 100. The fact that there is 
no simple scaling relation with /nl supports our concerns regarding 
the separability of lensing induced and primordial non-Gaussianity. 
To some extent the /NL-dependence of the bias qualitatively ap- 
proaches a quadratic form in the acoustic regime and beyond. 

It is important to keep in mind that these are just qualitative 
considerations. The actual discrepancy between the bias obtained 
from a rigor simultaneous treatment of primordial as well as lensing 




Figure 3. Additional bias associated with the presence of primordial non- 
Gaussian curvature perturbations. The bias (red curve) has been extracted 
from numerical simulations according to equation (37) with /nl = 100. In 
addition, the analytical estimates of the bias (green and blue curve), phe- 
nomenologically amplified, are shown. 




Figure 4. Same as Figure 3 but for another choice of the /nl parameter, 
namely /nl = 1000. Note that in this plot the term has been scaled 
differently highlighting that the bias found in numerical simulations does 
not exhibit a simple scaling relation with /nl in contrast to the analytical 
model. 



induced non-Gaussianity via numerical simulations and its analyti- 
cal lowest order estimate is tremendous. The analytical expression, 
even phenomenological amplified, is not suited for an accurate and 
consistent bias subtraction scheme. However, and this is the most 
important result, for intermediate values of / N l, as suggested by 
current observational data, the bias associated with primordial non- 
Gaussianity of local type does not effect CMB lensing power re- 
construction at all. 



5 SUMMARY 

In this paper we revisited the influence of primordial non- 
Gaussianity on CMB lensing power reconstruction. We determined 
the bias associated with the non-trivial structure of the CMB 
trispectrum generated by non-Gaussian curvature perturbations of 
local-type. While previous work of Lesgourgues et al. (2005) ad- 
dressed this issue analytically, we used numerical simulations. This 
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allows to overcome the non-obvious hypothesis that both effects 
leading to a connected part of the CMB four-point function do 
completely separate. The results of this work can be summarized 
as follows: 

(i) The analytical expressions derived by Lesgourgues et al. 
(2005) in the limit of a flat sky have been extended to the all-sky 
formalism. 

(ii) In the case of non-Gaussian primordial curvature perturba- 
tions of local type the bias arising from a non-trivial trispectrum of 
the intrinsic CMB temperature contrast is substantially larger than 
so far expected. 

(iii) The limitation of the simulated non-Gaussian CMB maps 
(discussed in Section 4.2) does not change the conclusion. Discard- 
ing the physical mechanism and just thinking of the simulated maps 
as the simplest example of a scale dependent form of non-Gaussian 
structure in the unlensed CMB, one finds that the associated con- 
nected part of the unlensed CMB is potentially a severe contami- 
nant in lensing power reconstruction studies. This interpretation of 
the simulated data addresses another important aspect. In order to 
estimate the amplitude of the bias prior knowledge of the type of 
primordial non-Gaussianity is needed. For removing the bias this 
prior knowledge is even more important. 

(iv) The strength of the contamination is obviously set by the 
amplitude of the non-Gaussian contributions. For current values of 
/nl only smallest multipoles are affected, while the effect on in- 
termediate and small scales is completely negligible. Low multi- 
poles are anyway cosmic variance dominated and have therefore 
only small predictive power. Thus, even in the full-sky case the 
bias arising from primordial non-Gaussianity of local type is not 
expected to compromise CMB lensing power reconstruction at all. 

As this work focuses on the bias associated with primordial non- 
Gaussian curvature perturbations solely characterized by /nl an 
extension of this study including also the contributions from the 
gNL-t erm is desirable. However, the simulation of accurate CMB 
maps containing the correct non-Gaussian structure arising from 
both, /nl and gNL, is challenging. 
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APPENDIX A: SOME FLAT-SKY EXPRESSIONS 

In this appendix we compile the flat-sky expressions used in this 
work. Starting with the lensing reconstruction bias linear in the 
lensing potential (cf. equation 12), the normalization reads 

A-\L) = J ^iL f(\ u \ 2 )g{\ u \ 2 ), (Al) 

while the optimal weights are given by 

/( 1„1 2) l.-LC^ + VLCfQ 
8W' l 2> - ?r00 r @@ ~ 7r oo r ee {A2 > 

"^A.obs^.obs ^Wi.obs^.obs 

with the additional condition li + I2 = L imposed. 

Finally, the expression for the intermediate quantities 'Pj'^(L) 
which build up the flat-sky fully reduced CMB trispectrum associ- 
ated with primordial non-Gaussianity of local type (cf. equation 24) 
reads 

= J r]dr { J r 2 2 dr 2 F L ( ri ,r 2 ) 



x (a e3 (r 2 )j3 ei (r 2 ) + a (i (r 2 )/3 { ,(r 2 )) (A3) 
with F L (fi, r 2 ), a t (r) and p[(r) defined in equations (17) - (19). 



This paper has been typeset from a TjX/ ETpX file prepared by the 
author. 



